EURONU-WP6-09-12 
IDS-NF-010 

Neutrino factory in stages: Low energy, high energy, off-axis 



Jian Tang & , Walter Winter 13 

a ' b Institut fur Theoretische Physik und Astrophysik, Universitat Wiirzburg, 

D-97074 Wiirzburg, Germany 



Abstract 

We discuss neutrino oscillation physics with a neutrino factory in stages, including the possi- 
bility of upgrading the muon energy within the same program. We point out that a detector 
designed for the low energy neutrino factory may be used off-axis in a high energy neutrino 
factory beam. We include the re-optimization of the experiment depending on the value of 
#13 found. As upgrade options, we consider muon energy, additional baselines, a detector 
mass upgrade, an off-axis detector, and the platinum (muon to electron neutrino) channels. 
In addition, we test the impact of Daya Bay data on the optimization. We find that for large 
$13 ($13 discovered by the next generation of experiments), a low energy neutrino factory 
might be the most plausible minimal version to test the unknown parameters. However, if 
a higher muon energy is needed for new physics searches, a high energy version including 
an off-axis detector may be an interesting alternative. For small #13 (#13 not discovered by 
the next generation), a plausible program could start with a low energy neutrino factory, 
followed by energy upgrade, and then baseline or detector mass upgrade, depending on the 
outcome of the earlier phases. 
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1 Introduction 



Three-flavor neutrino oscillations have been accepted as the most successful interpretation of 
neutrino flavor changes, see, e.g., Ref. [1]. In particular, the solar and atmospheric oscillation 
parameters have been measured with very high precisions, and the reactor mixing angle 9\ 3 
has been strongly constrained. Future long-baseline and reactor neutrino experiments will 
test this small angle further, and be sensitive to leptonic CP violation and the neutrino mass 
hierarchy (see Ref. [2] and references therein). The ultimate instrument for these purposes 
might be a neutrino factory [3-6]. Using different baselines and oscillation channels, it 
can basically disentangle all of the remaining oscillation parameters [7-10] in spite of the 
presence of intrinsic correlations and degeneracies [6, 11-13] for extremely small values of 
#13. Furthermore, a neutrino factory and other future neutrino oscillation experiments will 
be sensitive to new physics searches, see Ref. [2] for a summary. 

The design of a neutrino factory has been put forward and discussed in international studies, 
such as in Refs. [2,14-16]. Especially the most recent study, the International Neutrino 
Factory and Superbeam Scoping Study [2, 17,18], has laid the foundations for the currently 
ongoing Design Study for the Neutrino Factory (IDS-NF) [19]. This initiative from about 
2007 to 2012 is aiming to present a design report, schedule, cost estimate, and risk assessment 
for a neutrino factory. It defines a baseline setup of a high energy neutrino factory (HENF) 
with = 25 GeV and two baselines L x ~ 4 000 km and L 2 — 7 500 km (the "magic" 
baseline) operated by two racetrack-shaped storage rings, where the muon energy is 25 GeV 
(for optimization questions, see Refs. [5,6,10,20-24]). A key component is the magnetized 
iron detector (MIND) as far detector, where the magnetization is necessary to distinguish 
the "right-sign" (e.g., from — > v^) from the "wrong-sign" (e.g., from v e — > 9^) muons. 
The identification of the muon charge of the wrong-sign muons allows, for example, for 
CP violation measurements in the muon neutrino appearance channels [6,20]. Possible 
near detector configurations for cross section and flux measurements have been discussed in 
Refs. [17,25]. 

As a more recent development, a low energy neutrino factory (LENF) with ~ 4 GeV to 
5 GeV has been proposed as an alternative to the HENF [26-29] . The main purpose of this 
alternative has been the reduction of accelerator cost in the case of large #13. While the 
high energy neutrino factory relies on the MIND, the low energy neutrino factory as been 
proposed with a magnetized Totally Active Scintillator Detector (TASD), which allows for 
a lower threshold, better energy resolution, and (possibly) electron charge identification, 
which is required for the so-called "platinum" (v^ — > v e and — > v e ) channels. These 
channels are the T-inverted channels of the muon neutrino appearance channels. Because 
of the same matter effects in the — > u e and v e — > (or — > V e and u e — > u^) channels, 
CP violation can, in principle, be extracted without convolution with the matter effects. 1 
For the LENF, the useful number of muon decays may be increased by about 40% by an 
optimization of the neutrino factory frontend for low energies [29]. Since the spectrum of 
a LENF can also be produced by a HENF in an off-axis detector (OAD), we consider this 

1 In the CP-conjugate channel, the matter effects arc different, which means that the fundamental CP 
violation has to be disentangled from the Earth matter effects, which violate CP extrinsically (Earth matter 
does not contain any antimatter); see, e.g., Ref. [30]. 
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option as well. Such a detector would have a significant increase in low energy events. 

As far as possible physics outcomes are concerned before the decision for a neutrino factory, 
we restrict ourselves to standard oscillation physics. The small mixing angle #13 may be 
discovered by upcoming reactor (such as Double Chooz and Daya Bay) and superbeam (such 
as T2K and NOz/A) experiments until about 2012 to 2015 if sin 2 26 13 > 0.01 [31,32], whereas 
for sin 2 29\ 3 < 0.01, there will only be a new exclusion limit. At around this time, a reference 
design report will be available by the IDS-NF [19], possibly with two different setups (a 
LENF and a HENF). In addition, LHC may have sufficient luminosity to indicate the new 
physics scale, which may point towards a muon collider, in favor of a neutrino factory. The 
decision for any future facility has to be based on this knowledge (plus additional potential 
new physics discovered in oscillation and non-oscillation searches). We therefore refer to 
the sin 2 29i 3 ~ 0.01 as the splitting point between the "small #13" and "large $13" cases, 
which marks the end of the reactor and (first generation) superbeam experiment dominated 
era. We use the terms u 9 13 not found by the next generation" and "6*i 3 found by the next 
generation" equivalently. For the small 9i 3 case, the region 0.01 < sin 2 2#i 3 < 0.1 therefore 
is referred to as "next generation excluded". In our work, the large 6>i 3 case is treated in 
Sec. 3, the small #13 case in Sec. 4. 

For several reasons, such as external boundary conditions, a neutrino factory complex may 
not be built at once, but instead be regarded as step-by-step program, perhaps, towards a 
muon collider. In this study, we discuss both the LENF and HENF, where we are particularly 
interested in upgrade scenarios. This means that we demonstrate how building a neutrino 
factory in stages makes sense physics-wise, and we illustrate how the knowledge from earlier 
data affects the optimization. For the large #13 case (Sec. 3), we discuss the minimal 
requirements for a neutrino factory setup to measure the yet unknown parameters for both 
LENF and HENF. In addition, we re-consider the baseline optimization for different matter 
density uncertainties, the presence of the platinum channel or an off-axis detector, or the 
inclusion of Daya Bay data, which will be dominating the sensitivity to #13 at that time [32] . 
Furthermore, we show how much adding a second baseline or the combination between 
LENF and HENF will buy. For the small #13 case (Sec. 4), statistics in an off-axis detector 
will be very small, and platinum (because of charge identification backgrounds) or Daya 
Bay data (because of statistics) will not contribute significantly. In this case, we show 
a conceptually plausible upgrade scenario starting with a LENF, followed by an energy 
upgrade, and then by another detector upgrade or second baseline. We illustrate how the 
preceding phases affect the optimization, and we discuss possible synergies between LENF 
and HENF. In the following Sec. 2, we first of all motivate our simplifying assumptions, and 
we describe the experiments and possible upgrades. 

2 Assumptions and description of experiments and upgrades 

Here we describe our simplifying assumptions, experiments used, possible upgrades we will 
discuss, and the simulation techniques. 
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2.1 Simplifying assumptions 



For the sake of simplicity, we make the following assumptions for this study: 

• We only discuss two different muon energies = 4.12 [27] (LENF) and = 
25 GeV [19] (HENF). The baseline-muon energy optimization has been performed 
in Refs. [10,28]. For all performance indicators, the conclusion has been that the 
desired performance can be reached with a certain threshold muon energy, but higher 
muon energies do not harm (within reasonable ranges). The muon energy choice of the 
HENF was driven by this observation, and = 25 GeV was as reasonable minimum 
for the MIND detector. In addition, note that the sensitivity to non-standard effects 
saturates at about this energy, because the matter resonance in the Earth's mantle 
can be covered at the peak [24]. The LENF choice was motivated by the large #i 3 
performance given the TASD detector, with the boundary condition of low accelera- 
tion effort. More recently, somewhat larger E^ ~ 4.5 [29] are considered, which do 
not affect this discussion qualitatively. 

• We use the TASD as off-axis detector in the HENF with = 25 GeV with an off- 
axis angle of 0.55°. The considered spectrum of this off-axis neutrino factory detector 
(OAD) is practically identical to that of the LENF for this off-axis angle, which is a 
feature of the neutrino factory flux (see, e.g., Ref. [25]). Since the magnetized TASD 
is the detector used for the LENF and (so far) only discussed for that beam spectrum, 
we use the same beam spectrum. For E^ = 25 GeV, the corresponding off-axis angle 
is about 0.55°, or about 38 km for L = 4 000 km. In summary, we only use the MIND 
detector for HENF (on- axis), and the TASD for LENF (on- axis) and HENF (off-axis). 

• We normalize the luminosity of the neutrino beam produced by one decay straight 
to that of the IDS-NF baseline setup. The IDS-NF setup uses 2.5 ■ 10 20 useful muon 
decays per year, polarity, and decay straight, which corresponds to scale factor (SF) 
one. Since this setup uses two storage rings, the muons have to be shared among the 
rings. Therefore, if only one storage ring is used, such as for a one-baseline LENF or 
HENF, or a HENF+OAD, all muons can be injected in the same ring, and we have 
SF=2. If two baselines are operated simultaneously, on the other hand, we have SF=1. 
Note that the currently discussed LENF uses SF=2.8 [29], where the additional 40% 
increase comes from an optimization of the frontend. 2 Therefore, we sometimes show 
this as nominal scale factor. However, our approach has the advantage that LENF 
and HENF-OAD have the same performance. In addition, we require that the same 
frontend be used for LENF and HENF, which seems to make similar luminosities 
plausible. 

• We assume that the platinum channel is only relevant for the TASD, i.e., LENF or 
HENF-OAD. The reason is that the electrons create showers, which leads to an upper 
threshold for platinum charge identification (and possibly unpredictable behavior). 

2 Thcse numbers are based on 10 7 s operation time of the accelerator at the nominal luminosity per year. 
For 2 • 10 7 s operation, one has SF=5.6 for the LENF, as the more aggressive setup in Ref. [29]. We use 
10 7 s for all experiments in this study. 
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• We focus on standard three flavor oscillations with the three standard performance 
indicators: #i 3 discovery, mass hierarchy (MH) discovery, and CP violation (CPV) 
discovery. 

2.2 Low energy neutrino factory and possible upgrades 

Our LENF is based upon Ref. [27]. It uses = 4.12 GeV and a magnetized Totally Ac- 
tive Scintillator Detector (TASD) with 20 kt fiducial mass times efficiency at (typically) 
L = 1300 km. We will re-consider the optimization of the baseline including different types 
of upgrades. The detector threshold is assumed to be 500 MeV and the detection efficiency 
73%. Compared to Ref. [27], we use an energy resolution AE — 0.1 ^/E/GeV GeV, which 
is typical for TASD detectors such as NOz/A. Both muon neutrino (and antineutrino) dis- 
appearance and appearance channels are included. The background level is conservatively 
estimated to be 10~ 3 . Note that we include two types of backgrounds for the appear- 
ance channels, one which scales with the disappearance rates (such as from charge mis- 
identification), and one which scales with the un-oscillated spectrum (such as from neutral 
current events), both at the level of 10~ 3 . As in Ref. [27], we take the systematical errors 
to be 2% for all signal and background errors. 

As potential "upgrades" , these may the most relevant options: 

Platinum channel Our description of the platinum — > v e and — > u e channels is based 
upon Refs. [10,33,34]. The efficiency is 40%, the energy resolution is AE = 0.15 E, and 
the background level is about 0.01. Note that the platinum channel is very difficult 
to use for higher energies, since the electrons induce showers for which the charge 
is difficult to measure. Therefore, we only use it in the TASD for the LENF (and 
HENF-OAD). Because of the relatively high background level, we only consider the 
platinum channel for large 9 13 (c/., Ref. [10]). 

Energy upgrade We consider an energy upgrade to the HENF (see below). Note that 
since the higher muon energies require a different detector technology, we use the 
MIND detector in that case. 

External 6 13 measurement For large 9 13 , we consider the impact of an external mea- 
surement of 13 . Since we expect that Daya Bay dominates the #i 3 sensitivity at the 
time given, we only use the data from this experiment simulated as in Ref. [32]. 

During the completion of this work, Ref. [29] has appeared, which uses = 4.5 GeV and 
energy-dependent detection efficiencies. The background level for muon neutrino appearance 
and platinum channels in Ref. [29] is effectively a factor of two smaller than in this study. 
The energy resolution is similar (10%, linear in E). We have checked that the differences in 
efficiencies and energy resolution have practically no impact on the sensitivities. However, 
we conservatively assume higher background levels and a lower luminosity, which do have 
some impact. Therefore, our LENF has to be understood as conservative "minimal" choice 
LENF. As we will demonstrate, this choice is sufficient for the measurements for large 
sin 2 20i3 • 
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2.3 High energy neutrino factory and possible upgrades 

Our HENF is based upon the current IDS-NF baseline setup. In the standard configura- 
tion, it uses two 50 kt magnetized iron detectors (MIND) at baselines of about 4000 km and 
7500 km. We start with one baseline first (4000 km for small #13), and then consider poten- 
tial upgrades. Both disappearance and appearance channels are included. The energy 
resolution of the detectors is AE = 0.55^E/GeV GeV. The disappearance channels 
have a 90% detection efficiency with a threshold of 1 GeV. The detector efficiencies (up to 
74%) and background fractions (between about 10~ 5 and 10~ 4 ) of the appearance channels 
are a function of the neutrino energy; for details, see Ref. [19]. The systematical errors for 
signal and background normalizations are taken to be 2.5%, uncorrelated among all chan- 
nels. The number of useful muon decays per year, polarity, and storage ring is 2.5 10 20 (for 
two baselines operated simultaneously), which corresponds to SF=1. The simulation of the 
neutrino factory is based on Refs. [10,13,19,35]. 

As potential "upgrades" , we consider the following options: 

Off-axis detector A potential OAD is considered at the off-axis angle 0.55°, where the 
beam spectrum corresponds to the LENF spectrum. Therefore, the TASD from the 
LENF is used at the off-axis site. Of course, the platinum channel may be used in 
this detector as an additional upgrade. 

Second baseline A second baseline is already included as upgrade in the IDS-NF baseline 
setup. Of course, it requires an additional MIND, which we assume to have 50 kt. 
As baseline, typically 7500 km is used for small 6>i 3 . However, the HENF (and also 
LENF) optimization for large #13 might be different, and will be studied. Note that 
if both baselines are operated simultaneously, fewer muons can be injected in each 
storage ring. 

Detector mass upgrade For small $13, a detector mass upgrade of the MIND detector 
from 50 kt to 100 kt at the shorter baseline will be considered. Such an upgrade is 
currently discussed already within the IDS-NF. For large $13, an upgrade of the TASD 
would disproportionally increase the detector cost compared to the whole complex. 

Knowledge of matter density profile For large #13, the matter density uncertainty is 
known to be one of the limiting factors (see, e.g., Refs. [10,13,36]), together with sys- 
tematics. Therefore, we test the impact of the matter density knowledge. This implies 
that we also test how much the performance would improve (or how the optimization 
would change) in the presence of different matter density uncertainties. 

2.4 Simulation techniques 

The total running time of our experiments is assumed to be ten years, unless stated oth- 
erwise. For the useful number of muon decays, we use the "scale factor" SF, where SF=1 
corresponds to 2.5 • 10 20 useful muon decays per year, which is the IDS-NF standard per 
baseline and polarity. In one storage ring, both polarities (/i + and /i - decays) are operated 
simultaneously. If only one baseline is operated, SF=2, because all muons are injected in 
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the same storage ring. The scale factor re-scales the total luminosity, and can be used to 
test the impact of any re-scaling of statistics. Therefore, we sometimes use it as parameter. 

For the sensitivity analyses we use the oscillation parameter values (see, e.g., Refs. [1, 
37]): Am 2 21 = 7.65 • 10" 5 eV 2 , \Am 2 31 \ = 2.40 • 10" 5 eV 2 , sin 2 12 = 0.304, sin 2 6 23 = 0.500, 
and a normal hierarchy, unless stated otherwise. We impose external la errors on Am| x 
(4%) and 0u (4%) as conservative estimates for the current measurement errors [37]. We 
do not include an external measurement of the atmospheric parameters. In addition, we 
include a 2% matter density uncertainty, unless stated otherwise [36,38]. For the experiment 
simulation, we use the GLoBES software [39,40]. All "unused" oscillation parameters are 
marginalized over, such as all parameters except for sin 2 26>i 3 for the sin 2 2#i 3 discovery 
reach. 



3 Neutrino factory for large #i 3 (#i 3 found) 

In this section, we discuss the strategy for a neutrino factory if sin 2 2#i 3 is found by the next 
generation of experiments, i.e., sin 2 2#i 3 > 0.01. We first of all discuss the minimal require- 
ments for a neutrino factory in terms of baseline and luminosity without any additional 
information or upgrades. Then we study the baseline optimization including upgrades. 
Finally, we discuss two-baseline configurations for large #i 3 . 

3.1 Minimal neutrino factory 

Here we investigate the minimal requirements for a neutrino factory. We consider both the 
LENF and HENF, but keep in mind that "minimal" clearly refers to the LENF. In addition, 
we only consider one baseline. We assume that the measurement should be independent of 
that from other experiments, such as Daya Bay, and we do not consider any additions or 
upgrades. 

Compared to the small sin 2 2# 13 case, it is much easier to define a minimum wish list if 
sin 2 2# 13 has already been observed. Here we follow the minimum wish list in Ref. [41] 
(which was discussed there in context of the beta beam): 

1. 5cr independent confirmation of sin 2 2# 13 > (for any Sep). 

2. 3<t determination of the mass hierarchy (MH) for any (true) Sep. 

3. 3(7 establishment of CP violation (CPV) for a certain fraction (such as 80%) of all 
(true) 5 C p. 

The only "arbitrary" in this list is the fraction of 5 C p for which CPV should be discovered. 
A fraction of 80% corresponds to Cabibbo-angle precision, which can be motivated in quark- 
lepton complementarity scenarios (see, e.g., Refs. [42,43]). Alternatively, it corresponds to 
the precision of the CP phase in the quark sector is measured. In this wish list, point 1 
is typically easy for most of the parameter space; therefore, we do not show it explicitely 
anymore. Point 2 typically requires a certain minimum baseline. Point 3 requires sufficient 
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Figure 1: Discovery of CPV (dark/red) and MH (medium gray/light blue) for the one baseline (minimal) 
LENF as a function of baseline and luminosity scale factor SF. Discovery reach is given within the shaded 
regions at the 3er CL, where for CPV a fraction of 5qp of 75% or 80% is required (as indicated), and for 
the MH a fraction of 5qp of 100%. The stars show the baseline with the minimal SF: in the left panel 
(1100 km, 2.0) and in the right panel (1150km, 2.6). The nominal luminosity is given by SF=2.8. Here the 
true value of sin 2 29\s is chosen as given in the plot panels, and a normal hierarchy is assumed. The matter 
density uncertainty is assumed to be 2%. 

luminosity and an appropriate baseline window close to the oscillation maximum. Compared 
to Ref. [41], we simplify the analysis somewhat and show the results only for particular 
choices of (true) sin 2 2^13 (not ranges allowed by the next generation of experiments). Of 
course, the choice of sin 2 29 13 will be motivated by the results from preceding experiments. 
In addition, we show the normal hierarchy only. 

In order to identify the minimal version of the neutrino factory, we re-optimize the baseline, 
and, at the same time, identify the minimum luminosity for the optimal baseline with respect 
to the above performance indicators. For the LENF, we show in Fig. 1 the discovery reach for 
CPV (dark/red) and MH (medium gray/light blue) as a function of baseline and luminosity 
scale factor SF. Discovery reach is given within the shaded regions at the 3a CL, where for 
CPV a fraction of 5qp of 75% or 80% is required (as indicated), and for the MH a fraction 
of <5cp of 100%. Sensitivity to both performance indicators is given in the overlap region, 
where one typically also has #13 discovery potential for all Sep- The stars mark the points 
with the minimal SF where all performance indicators can be measured; they therefore show 
the "minimal configurations". For the LENF, the minimal baseline is determined by the 
MH reach, and the minimal SF by the CPV reach. The nominal luminosity (SF=2.8) is 
sufficient for the CPV measurement for 80% of all true 5qp and for the MH measurement 
for all 6qp in the baseline window 1100 km < L < 1400 km for both values of sin 2 26* 13 
(left and right panel). One can read off these figures that luminosity is clearly an issue for 
large sin 2 29^. If, for instance, only a lower SF can be achieved, the CPV discovery reach 
decreases accordingly. The "minimal" (optimal) LENF, i.e., the one with the lowest SF, is 
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Figure 2: Discovery of CPV (dark/red) and MH (medium gray/light blue) for the one baseline (minimal) 
HENF as a function of baseline and luminosity scale factor SF. Discovery reach is given within the shaded 
regions at the 3cr CL, where for CPV a fraction of Sqp of 75% or 80% is required (as indicated), and 
for the MH a fraction of <5cp of 100%. The stars show the baseline with the minimal SF: in the left 
panel (3500 km, 1.35), in the middle panel (2100km, 2.8), in the right panel (5400 km, 4.8). The nominal 
luminosity is given by SF=2. Here the true value of sin 2 26 l 13 = 0.08 and a normal hierarchy are assumed, 
and the matter density uncertainty is varied from the left to the right, as shown in the plot panels. 

in both panels at about L ~ 1100 km. However, the FNAL-DUSEL baseline L = 1290 km 
is close enough to optimum. 

For the HENF, we show in Fig. 2 the discovery reach for CPV (dark/red) and MH (medium 
gray/light blue) for the one baseline (minimal) HENF as a function of baseline and luminos- 
ity scale factor SF. Here we show the dependence on the matter density uncertainty, because 
this is the main impact factor for the HENF, and instead choose only one value of (true) 
sin 2 26* 13 = 0.08. Compared to the LENF, the CPV discovery reach typically determines the 
minimal configuration denoted by the stars. One can easily see that both performance and 
baseline optimization strongly depend on the degree the matter density profile is known. 
For instance, in the left panel, Ap/p = (perfectly known matter density profile), and a 
fraction of 5cp of 80% can be easily achieved in a baseline window 2000 km < L < 5000 km 
for the nominal SF=2, where also the MH can be measured for all 5cp- For more real- 
istic matter density uncertainties, however, the situation becomes more complicated (c/., 
middle and right panels for 2% and 5%, respectively). The CPV discovery reach, in these 
cases, exactly deteriorates in the window which is optimal for no matter density uncer- 
tainty, and two local minima, one at a shorter baseline L ~ 2000 km and one at a longer 
baseline L ~ 5500km, remain (c/., "7r-transit" problem [13]). The absolute reach is in both 
cases above the threshold of 80%. Therefore, a reliable optimization of the HENF without 
upgrades for large sin 2 2#i3 is only possible if the matter density profile is precisely known. 

Since for realistic matter density uncertainties the LENF outperforms the HENF at the 
nominal luminosity and the HENF relies on a higher effort on the accelerator side, the LENF 
is probably the best "minimal" version of a neutrino factory for large sin 2 26*13. However, it 
may be upgraded for different purposes later. For instance, the search for new physics, such 
as non-standard interactions, requires higher muon energies [24]. Therefore, we nevertheless 
consider the HENF for large sin 2 2#i3, especially in the context of upgrades. 
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Figure 3: Fraction of Sep for which CPV will be discovered as a function of baseline L for the LENF 
(left panel) and HENF (right panel); 3cr CL. The different curves show the impact of the matter density 
uncertainty, as given in the plot legends. The vertical lines mark the regions where also the MH can be 
determined for all Sep at the 3cr CL. Here sin 2 2#i3 = 0.04, SF=2, and a normal hierarchy is assumed. 

3.2 Single baseline neutrino factory: Upgrades and optimization 

Here we discuss the LENF and HENF baseline optimization for large #13 together with 
possible upgrades or external input. In this case, we choose the nominal SF=2 for both 
LENF and HENF. This choice has the advantage that LENF can also be interpreted as 
HENF-OAD. 

As the first aspect, we show in Fig. 3 the fraction of Sqp for which CPV will be discovered 
as a function of baseline L for the LENF (left panel) and HENF (right panel). The different 
curves show the impact of the matter density uncertainty, as given in the plot legends. 
For the LENF (left panel), we find that the fraction of 5qp drops sharply off the optimal 
baseline window identified above. The matter density uncertainty hardly has any impact 
on the absolute performance or optimization. A fraction of 5qp of 80% can be reached 
at about the peak without any further upgrades. Note again that the LENF could also 
be HENF-OAD (therefore we chose SF=2 consistently), and option which we will discuss 
below. For the HENF (right panel), the optimal baseline depends on the matter density 
uncertainty, and also somewhat on sin 2 29^ (here we only show one example). The fraction 
of £cp of 80% can (for sin 2 29 13 = 0.04) be reached if the matter density is small enough 
(< 2%). For the chosen value of sin 2 2# 13 , the optimal baseline for both CPV and MH is 
around 2000km to 4000km (depending in matter density uncertainty and sin 2 2# 13 ). Again, 
the LENF is the more robust version for large sin 2 26*13, because the optimization does not 
depend on the matter density uncertainty. 

As far as possible upgrades are concerned, we show in Fig. 4 the fraction of Sqp for which 
CPV will be discovered as a function of baseline L for the LENF (left panel) and HENF 
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Figure 4: Fraction of Sep for which CPV will be discovered as a function of baseline L for the LENF (left 
panel) and HENF (right panel); 3er CL. The different curves show different upgrade options: off-axis TASD 
(OAD) and platinum channel. The thick curves give the baseline range where also the mass hierarchy can 
be determined for any Sep at the 3a confidence level. Here sin 2 29^ = 0.04, SF=2, a normal hierarchy and 
a 2% of the matter density uncertainty are assumed. 

(right panel), where we have chosen a matter density uncertainty of 2%. The thick curves 
give the baseline range where also the mass hierarchy can be determined for any 5qp at 
the 3cr confidence level. For the LENF (left panel), the most plausible upgrade may be 
the platinum channel, which clearly increases the absolute performance for large sin 2 26* 13 
(left panel) up to a fraction of <5cp 85%. In addition, the optimal baseline window becomes 
larger, and a short baseline, such as FNAL-Soudan, could be sufficient to break the mass 
hierarchy degeneracy (where the short baseline cutoff in the CPV discovery reach comes 
from) . 

A different possibility together with upgrades my be the HENF (right panel). An additional 
OAD would combine the virtues of the LENF and HENF. In this case, the optimal baseline 
becomes shorter, possibly as short as FNAL-DUSEL. Together with the platinum channel 
(in the OAD), CPV can be even measured for almost 90% of all Sqp at a baseline of about 
1100 km. At this baseline, the matter density uncertainty hardly has any impact. Note, 
however, that even longer baselines up to 4500 km still lead to a CPV discovery for more 
than 80% of all 5qp- In this case, the matter density profile has to be controlled at the level 
of a few percent. The correlation between the profiles seen by the off- and on-axis detectors 
(distance only about 38 km at L — 4000 km), which we have included in the figure, helps 
somewhat. Therefore, if for some reason (such as new physics searches) a HENF may be 
the preferred option, the combination with the OAD-TASD could be the most interesting 
upgrade. It allows for a relatively wide baseline window. 

We have also tested the impact of 5 year Daya Bay data, simulated in Ref. [32], which may 
affect the optimization for large #13. However, we have not found any significant impact, 
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Figure 5: Two baseline optimization for CPV: Fraction of Sep for which CPV is discovered as a function 
of the two neutrino factory baselines; 3a CL. Here the combinations LENF-LENF (left panel), HENF-LENF 
(middle panel), and HENF-HENF (right panel) are shown. Note that LENF can also refer to the OAD in 
this figure. Here a 2% matter density uncertainty, SF=1 (each baseline), the true value sin 2 26*i3 = 0.08, 
and a normal hierarchy are assumed. Here no upgrades are included. 

neither on absolute performance, nor on the baseline optimization, by the combination with 
the Daya Bay data. However, note that a certain simulated value of sin 2 29 13 is always 
assumed for the optimization, which is only known within the range expected from the next 
generation of experiments. 

3.3 Double baseline neutrino factory 

So far, we have assumed that only one baseline is used. A double baseline neutrino factory 
could also be an option for large #13, where the optimization might be different from the 
small #i3 case studied in the literature [9,24]. A double baseline neutrino factory could be 
interesting in terms of the following upgrades: 

Energy upgrade: LENF, followed by a HENF at a (possibly) different baseline 

Baseline upgrade: LENF at two different baselines or HENF at two different baselines 

Additional off-axis detector: HENF, followed by an OAD at a (possibly) different base- 
line 

We study these three options in Fig. 5 in terms of the fraction of 5cp for which CPV is 
discovered as a function of two baselines. For the baseline combinations, we have LENF- 
LENF, HENF-LENF (or OAD), HENF-HENF, where we choose SF=1 for each baseline. 
This means that we assume that two storage rings are required, and the muons are equally 
shared among the rings. In principle, LENF-LENF or HENF-HENF for the same baseline 
(along the diagonal) correspond to the previously studied cases. However, note that the 
external input (such as the knowledge on the solar parameters) is added twice here, and 
that the matter density uncertainty is assumed to be uncorrelated between the two baselines. 
Therefore, they cannot be directly compared to the previous results. 
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The main conclusion from Fig. 5 is that in all cases the LENF and HENF baseline optima 
are obtained almost uncorrelated with each other. This means that, if one combines two 
baselines and allows for two different storage rings, the LENF is optimal at about L = 
1000 km, the HENF at about L = 2000 km to 4000 km, which roughly corresponds to the 
individual one baseline optimizations. 3 Therefore, for any energy or baseline upgrade, one 
would follow the strategy from the individual baseline optimizations. For the combinations 
LENF+LENF and HENF+HENF, the most plausible upgrade may then be a detector mass 
upgrade. 

In summary, a LENF at a baseline of about 1100 to 1400 km may be the most plausible 
neutrino factory option for large sin 2 29 13. If possible, it should rely on electron neutrino 
appearance (platinum channel) as well. We pointed out that if new physics searches indicate 
that a HENF may be desirable, the LENF detector could also be used as OAD for the HENF, 
with significant impact on the sensitivity. 

4 Staging for small #i 3 (#i 3 not found) 

After several years of data taking from Daya Bay and the other next generation experiments, 
we will know whether sin 2 29 13 < 0.01 [32]. For this part, we therefore assume that 9 13 has 
not been found by these experiments, which corresponds to sin 2 2# 13 < 0.01. In this case, 
the combination with Daya Bay will not help very much. The same applies to the platinum 
channel, which is limited by the charge identification capabilites. How can one then build a 
neutrino factory step by step, while taking into account the knowledge from the preceding 
phases? 

What we already know is that L ~ 4000 km is optimal for CPV [10], L ~ 7500 km for a 
risk-minimized mass hierarchy measurement [10] and an excellent 6*13 sensitivity [9] (which 
approximately corresponds to the worst case discovery reach). The combination is good 
for all performance indicators [24]. However: These are optimized in the #13 direction, i.e., 
discovery for as small as possible #13. If we know, at some point, how large 6>i 3 is, the 
discussion may change, depending on the priorities we have. 

In this section, we therefore sketch a plausible staging scenario, starting with a low energy 
neutrino factory. As upgrade options, we consider increasing the muon energy, adding 
another baseline, and increasing the detector mass. For the sake of simplicity, we assume 
that the accelerator complex delivers 10 21 useful muon decays for both polarities per year. 
If there are two storage rings operated with the same muon energy, these muons can be (at 
least in principle) split arbitrarily between the rings. We proceed in three phases of data 
taking, five years each. Of course, there may be idle times between these phases, such as 
because of construction. Phase I represents a low energy neutrino factory, phase II includes 
the energy upgrade to a high energy neutrino factory, and phase III considers additional 
upgrades, such as a larger detector or an additional baseline. In any phase, we combine the 
data with the preceding phase. 

3 Only in some corners of the parameter space, such as HENF+HENF for sin 2 29\z = 0.12 and 5% matter 
density uncertainty, we have found that combining a shorter with a longer baseline may somewhat help. 
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Figure 6: Phase I: 9 discovery reach (3cr) as a function of baseline (left panel) and all discovery reaches 
(3a) as a function of sin 2 2#i3 and fraction of <5cp for the optimal 613 baseline L — 900 km (right panel). In 
the left panel, the best case (lower curve), worst case (upper curve) and "typical" (thick curve) true Sqp is 
shown. The "typical" <5cp corresponds to the median, i.e., the performance is better for 50% of all <5cp and 
worse for 50%. Here a normal hierarchy is assumed. 

4.1 Phase I: Low energy neutrino factory 

If #i3 is not discovered by the next generation of experiments, we assume that in phase I 
of a neutrino factory a low energy version with = 4. 12 GeV and a magnetized TASD 
as detector is operated. The main priority will be the search for # 13 , which means that 
the machine should be optimized for that. Our phase I has five years of operation with 
10 21 useful muon decays per year in both polarities, i.e., SF=2. We show in Fig. 6, left 
panel, the baseline optimization for the #13 discovery reach. In this figure, the best case 
(lower curve), worst case (upper curve) and "typical" (thick curve) true 6cp is shown. The 
"typical" 5qp corresponds to the median, i.e., the performance is better for 50% of all 
Sep and worse for 50%. Obviously, a baseline of about 900 km is close-to-optimal for the 
typical 5qp (thick curve), as it was for for CP violation for large (9 13 in Ref. [28]. For a 
risk-minimized performance (upper curve), somewhat longer baselines are preferred. In 
summary, baselines between about 500 km and 1500 km are sufficiently good. We use, in 
the following, L = 900 km, which will allow for a sin 2 2#i3 discovery for values between 0.003 
and 0.0003 for the normal hierarchy, depending on the true value of Sqp- 

In Fig. 6, right panel, we show the mass hierarchy, CP violation and #13 discovery reaches of 
phase I as a function of # 13 and 5cp- The region, which is accessible to the next generation 
of experiments and covered by the previous section, is shaded in gray. One can easily read 
off this figure that CP violation is only accessible for a small fraction of the parameter space 
if 9 is is discovered in phase I, and the mass hierarchy can practically not be determined for 
sin 2 2^ 13 < 0.01. Therefore, even in the case of a 9 13 discovery, phase II will be most likely 
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Figure 7: Results from phases I+II. The different regions show the #13, MH, and CPV discovery reaches 
(3<t) as a function of sin 2 2^13 and fraction of <5cp- The dashed curves represent phase II alone. Here a 
normal hierarchy is assumed. 

needed. In case of a sin 2 2#i3 discovery, we list some obtainable precisions for sin 2 29i 3 for 
different values here (90% CL, computed for the true Sep = 7r/2): 

For sin 2 26 13 = 0.001 : 0.00048 < sin 2 29 13 < 0.0037 (1) 
For sin 2 26 13 = 0.005 : 0.0037 < sin 2 26 13 < 0.0072 (2) 
For sin 2 26 13 = 0.01 : 0.0081 < sin 2 29 13 < 0.013 (3) 

4.2 Phase II: Energy upgrade? 

If #i3 is not discovered in phase I, the most plausible upgrade might be an energy upgrade, 
which comes together with a new detector and a new baseline (including storage ring) 
for optimal sensitivities. 4 For the energy, we use the IDS-NF standard muon energy = 
25 GeV. It is typically discussed together with the MIND detector. We assume that phase II 
corresponds to five years of data taking again, and we combine the data with phase I. Note 
that all useful muons in phase II are used for the new baseline (SF=2), i.e., 10 21 useful 
muon decays per year in both polarities. 



4.2.1 #i3 not discovered in phase I 

The optimization for small #13 has been extensively studied in the literature, see e.g. Ref [10]. 
We follow the conventional strategy to aim for discovery reaches for as small as possible #13: 

4 See, e.g., Ref. [10]: At least for small 6*13, L = 900km is far away from optimal for the higher muon 
energy. Therefore, we do not consider placing the new detector at the L = 900 km baseline. 
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Figure 8: Baseline optimization for the phase II baseline for the combination of phase I and phase II. 
Here the fraction of (true) 5qp for CPV discovery (3<r) is optimized for. The blue (dark) curves are drawn 
for a possible LENF in phase II, the red (light) curves for a HENF (energy upgrade) in phase II. The dashed 
curves do not include degeneracies. The thick curves show the baseline ranges where also the mass hierarchy 
can be determined for any Sep at the 3tr confidence level. The baseline for phase I and the baseline range 
of the IDS-NF (for the shorter) baseline are marked. The different panels correspond to different values of 
sin 2 2#i3, as they may be discovered in phase I. Normal hierarchy assumed. 

L ~ 4000 km is a good choice for CPV and also has an excellent #13 discovery reach for 
most values of Sep- 

We show in Fig. 7 the result for this upgrade. With the upgrade, 6*13 will be discovered for 
sin 2 2#i3 as small as 10 -4 for any value of 5qp- In the case of a discovery, CPV and the mass 
hierarchy can be determined for a large fraction of the parameter space for sin 2 29\3 > 10~ 4 , 
beyond the next generation of experiments. If the MH or CPV are not found, we discuss in 
Sec. 4.3.1 further possible upgrades. 

Very interestingly, the combination with phase I helps to resolve the degeneracies. This is 
illustrated by the dashed curves, which are computed for phase II alone. Without phase I, 
the degenerate solutions severely affect the CPV and MH discovery reaches, whereas the 
phase I has just enough statistics to partially resolve the degeneracies. In this case, the 
small matter effects in phase I (because of the shorter baseline and lower energy) become 
an advantage, because the mass hierarchy degeneracy has a very different location. In fact, 
we have also tested using double luminosity in phase II alone, to check if there is a real 
synergy beyond the addition of statistics (c/., Ref. [44] for a more detailed discussion). In 
fact, we could identify such a synergy for the MH and for CPV for a part of the parameter 
space (the larger sin 2 26*i 3 values). 

4.2.2 #i3 discovered in phase I 

In this case, we re-consider the baseline optimization, depending on the (9i 3 found in phase I. 
For the sake of completeness, we study the physics of a second baseline with and without 
an energy upgrade, i.e., LENF or HENF in phase II. The result for the optimization of the 
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second baseline is shown in Fig. 8 for different best-fit values of sin 2 2# 13 , as they might 
be obtained in phase I; c/., Eqs. (f) to (3). The plots show the fraction of (true) 5cp f° r 
CPV discovery (3a) as a function of baseline. The blue (dark) curves are drawn for the 
LENF, the red (light) curves for the HENF (in phase II). The dashed curves do not include 
degeneracies. The thick curves show the baseline range where also the mass hierarchy can be 
determined for any 5cp at the 3a confidence level. Note that the matter density uncertainty 
is assumed uncorrelated with the first baseline. 

For small values of sin 2 2# 13 at the lower end of the range reachable in phase I (left panel), 
the mass hierarchy cannot be measured for all Sep with the LENF in phase II, and the energy 
upgrade performs much better for CPV. The baseline choice L 2 ~ 4000 km is precisely at 
the optimum for the HENF. Therefore, the strategy will be the same in the case of sin 2 2#i3 
not discovered in phase I: energy upgrade with a MIND at L 2 — 4000 km. Note that the 
fraction of 5cp for LENF peaks at about 900 km, as it was optimized before. 

For large values of sin 2 2# 13 at the upper end of the range reachable in phase I (right panel), 
an energy upgrade may not be necessary, but a significantly longer baseline for the LENF 
in phase II is preferable to determine the mass hierarchy: 1600 km < L 2 < 4000 km. The 
energy upgrade (HENF) buys about 10% in the fraction of 5 C p, where slightly shorter 
baselines than in the IDS-NF baseline are optimal: 2000 km < L 2 < 4000 km. 

For the intermediate case (middle panel), the CPV performance is significantly worse than 
for the high energy version. Therefore, an energy upgrade seems to be desirable, within 
about the same baseline range 2000 km < L 2 < 4000 km. Note that the LENF in phase II 
can measure the mass hierarchy for about 70% to 90% of all 5qp in the range 2000 km < 
L 2 < 4000 km, which is below threshold in the figure. 

Since L 2 ~ 4000 km is sufficiently close to optimum, we choose that baseline together with 
the energy upgrade for the following discussion in phase III, such as for CPV not found in 
phase II. 

4.3 Phase III: Third baseline or detector upgrade 

Here we consider two possible further upgrades: Another baseline for the HENF, or a 
detector mass upgrade from 50 kt to 100 kt for the MIND detector at the same location at 
L = 4000 km. If only one baseline is used, all muons are delivered to the respective storage 
ring. If a second baseline is added, the muons can be (almost) arbitrarily split between 
the two possibly simultaneously operating storage rings. We have tested a simultaneous 
operation of both baselines in phase III with an equal splitting of the muons between the 
two storage rings, and we have tested an option with all muons in the new storage ring. 
Since the differences are moderately small, we choose the second option. Again, the data 
from phases I and II are added. For the second baseline, we consider the magic baseline 
L = 7500 km, because it is known to lead to optimal performances for all performance 
indicators for small #13 in the combination with the shorter baseline [24] . 

We show the results for phase III in Fig. 9. In this figure, the solid curves correspond to the 
baseline upgrade in phase III (magic baseline), the dashed curves to a detector upgrade. 
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Figure 9: Results from phases 1+11+111. The different regions show the #13, MH, and CPV discovery 
reaches (3er) as a function of sin 2 26*13 and fraction of 5qp- The solid curves correspond to the baseline 
upgrade in phase III (magic baseline), the dashed curves to a detector upgrade. In addition, some error 
bars for #13 are shown for several selected best-fit points, as they are obtained from phase I+II. They are 
computed for the true Sqp — 7r/2 (as an example) at the 90% CL. Here a normal hierarchy is assumed. 

4.3.1 #i3 not discovered in phase I+II 

If #i3 is not discovered in phase I+II, the best option to increase the #13 is obviously in- 
creasing the detector mass in phase III (dashed curves in Fig. 9). In this case, The #13 and 
CPV reaches in the #13 direction can be significantly improved in phase III compared to 
phase I+II alone. That means that in the case of a #13 discovery, also the chance to observe 
CP violation increases. 

4.3.2 6 13 discovered in phase I+II 

If is discovered in phase I+II, we will know it with a certain precision before phase III. 
This is illustrated by the error bars in Fig. 9, which are shown for several selected best-fit 
points, as they could be obtained from phase I+II. Based on these results, the best upgrade 
strategy can be chosen. Since the MH and #13 can, in most cases, be easily measured, the 
only performance indicator to be optimized is the fraction of 5qp for which CP violation will 
be discovered. For instance, for sin 2 29 13 = 0.002 (best-fit), degeneracies are important, and 
the second baseline is the best strategy. This case also includes the possibility that sin 2 29\ 3 
is discovered already in phase I, but CPV has not been found in phase I+II. 

For sin 2 2^13 = 0.0001 (best-fit), however, the detector mass upgrade makes sense. For the 
other two cases shown in the figure, both options are equally good. Note that although 
the detector mass upgrade seems to be a good option in many cases, the resolution of 
degeneracies happens through sufficient statistics. This means that if the target luminosity 
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Figure 10: 6, MH, and CP discovery reaches for a neutrino factory in three phases: Phase I (light 
gray/green) is a low energy neutrino factory with a magnetized TASD, phase II (medium gray/red) adds 
an energy upgrade with a MIND at the 4000 km baseline, and phase III (dark gray/blue) includes another 
(magic) baseline (solid curves) or a detector upgrade at the 4000 km baseline (dashed curves) . 5a confidence 
level. 

cannot be reached or a higher confidence level is chosen, the magic baseline is the more 
robust solution. We illustrate this in Fig. 10 below for the 5er confidence level (where, for 
instance, the mass hierarchy discovery reach of the magic baseline option is superior). 

Note that, in principle, one could also optimize the second baseline again for an optimal 
CPV sensitivity depending on the #13 from phase I and II. However, here we focus on the 
magic baseline because it provides orthogonal, i.e., qualitatively different, physics, which 
may be also useful for non-standard measurements. 

We summarize the neutrino factory in stages for small sin 2 2^13 in Fig. 10 (5a), where a 
second baseline is used in phase III (solid curves) or the detector upgrade (dashed curves). 
At the 5cr confidence level, the strength of the magic baseline to resolve degeneracies becomes 
more pronounced, such as for the mass hierarchy discovery reach. In addition, there are 
many other applications of this baseline, see, e.g., Refs. [23,45-47], for non-standard physics, 
see, e.g., Refs. [24,48]. However, the largest increase in the discovery reaches will come from 
the energy upgrade. 

As far as the number of different storage rings is concerned, the LENF storage ring from 
phase I cannot be recycled in phase II since it points to a shorter baseline (for geometric 
reasons). Typically, it is also assumed to be much smaller because of the smaller gamma 
factor. In phase III a new storage is required if the magic baseline is chosen as option, 
whereas the detector upgrade does not require another storage ring. Therefore, either two 
or three different storage rings are needed in total. 
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5 Summary and conclusions 



We have discussed the optimization of a low energy neutrino factory (LENF) and high 
energy neutrino factory (HENF), including the possibility of an energy upgrade from LENF 
to HENF. We have also pointed out that the magnetized TASD, which is proposed in 
the context of the LENF, might be used as an off-axis detector (OAD) in the HENF, 
which exactly the same beam spectrum if the off-axis angle is chosen accordingly. We 
have tested the impact of luminosity, baseline length, platinum channel, Daya Bay data 
(as representative for the next generation of experiments), additional baselines, and the 
matter density uncertainty on the absolute performance and optimization. Conceptually, 
we have distinguished the small sin 2#i 3 case (sin 

2#i3 < 0.01) and the large sin 2 26 l3 case 
(sin 2 2^13 > 0.01), which correspond to #13 excluded or discovered by the next generation of 
experiments. 

For large sin 2 2# 13 , we have demonstrated that a moderate luminosity LENF at a baseline 
of about 1100 to 1400 km may be the most plausible "minimal effort" neutrino factory 
option to measure CPV and MH. Because of the short baseline and relatively low energies, 
it is robust with respect to the knowledge of the matter density profile, and it allows for 
additions or upgrades, such as electron neutrino appearance (platinum channel), which 
significantly increase the sensitivity. We have not found any significant impact on the 
performance if the data from a prior sin 2 2# 13 measurement, such as from Daya Bay, are 
directly added. However, the obtained rough value of sin 2 26 13 is helpful for the optimization 
(our optimization has been performed for certain true values of sin 2 2#i 3 ). 

We have also discussed the HENF for large sin 2 2#i 3 , which can basically perform the same 
measurements with a similar target sensitivity if the matter density can be controlled at 
a level below 2%. In this case, the optimal baseline would rather be 2000 km to 4000 km, 
depending on the value of sin 2 2# 13 found and the matter density uncertainty. An OAD 
detector, which may also allow for the platinum channel because of the lower energies in 
the off-axis spectrum, can significantly enhance the sensitivity. The reason is that the 
T-inverted golden channel corresponds to the platinum channel, which allows for a direct 
extraction of the intrinsic phase 5qp without extrinsic CP violation from the matter effect, 
compared to the CP-conjugated golden channel (such as the antineutrino channel). While 
the optimal baseline for CPV is at about 1000 km to 1500 km, baselines up to 4500 km still 
allow for a CPV measurement for 80% of all Sep- Note that a HENF could be interesting for 
large sin 2 29 13 for different reasons, such as a case for new physics searches at the neutrino 
factory. For example, for non-standard matter effects, high muon energies are mandatory 
for sensitivities beyond the current limits [24] . In this case, also the r production threshold 
is significantly exceeded, which allows for a detection of v T . For HENF+OAD, the baseline 
could be chosen in a wide window to search for new physics. 

For small sin 2 20 13 (sin 2 20 13 not discovered by the next generation of experiments), we 
have plotted a possible staging scenario starting with a LENF (phase I), followed by an 
energy upgrade (phase II) and then either a detector mass upgrade or second (magic) 
baseline (phase III). We have demonstrated that the LENF in phase I could significantly 
enhance the sin 2 2# 13 discovery reach beyond the exclusion limit expected from the next 
generation of experiments. The energy upgrade in phase II has the largest impact on the 
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sensitivities. Depending on the outcome of phase I (if sin 2 2# 13 is discovered), the baseline 
can be chosen accordingly, where larger values of sin 2 2# 13 prefer slightly shorter baselines 
than smaller values of sin 2 26 ) i 3 . We have also demonstrated that there is a synergy in 
resolving degeneracies between the LENF in phase I and the HENF in phase II. In phase III, 
one can either upgrade the detector mass, or the baseline. The decision for either of the 
two options can be based on the outcome of phase I+II. The advantage of the detector 
mass upgrade is that it does not require a new (large inclination) storage ring, and a single- 
baseline neutrino factory can be effectively operated with a factor of two higher luminosity 
than a two-baseline neutrino factory, because the muons do not have to be split between 
the different storage rings. The benefit of the magic baseline option is qualitatively different 
physics, such as relevant for new physics searches in the presence of more than six oscillation 
parameters (see, e.g., Ref. [48]). We summarize the neutrino factory in stages for small 
sin 2 2^ 13 in Fig. 10 (5a), where a second baseline is used in phase III (solid dark gray/blue 
curves) or a detector upgrade (dashed dark gray/blue curves). 

In summary, we have studied the optimization of a neutrino factory, where we have removed 
the constrained that the muon energy has to be fixed. We have found that for large sin 2 2#i 3 , 
a neutrino factory including a TASD (either on- or off-axis) has very good performance. If 
high muon energies are needed, such as to exceed the r production threshold, a HENF 
should include an off- axis detector. For small sin 2 29 13 , we have demonstrated that there 
can be a reasonable staging scenario including the LENF as first option, followed by an 
energy upgrade. We conclude that distinguishing LENF and HENF as clearly separate 
options may not be close to reality. A realistic program may include components of both 
options, no matter if sin 2 26*i 3 is small or large. 
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